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Abstract 

The shifted Tietz-Wei (sTW) oscillator is as good as traditional Morse potential in 
simulating the atomic interaction in diatomic molecules. By using the Pekeris-type ap¬ 
proximation to deal with the centrifugal term, we obtain the bound-state solutions of 
the radial Schrodinger equation with this typical molecular model via the exact quan¬ 
tization rule (EQR). The energy spectrum for a set of diatomic molecules (NO (a 4 II,:), 

NO (B * 2 n r ), NO (L' 2 cp), NO (& 4 S~), ICl (WS+), I Cl (A 3 rii) and ICl (A' 3 II 2 ) for 
arbitrary values of n and i quantum numbers are obtained. For the sake of complete¬ 
ness, we study the corresponding wavefunctions using the formula method. 
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1 Introduction 

By employing the dissociation energy and the equilibrium bond length for a diatomic 
molecule as explicit parameters, Jia et al [T] generated improved expressions for some well- 
known potentials including Rosen-Morse, Manning-Rosen, Tietz and Frost-Musulin poten¬ 
tial energy functions. These authors found that the well-known Tietz potential function is 
conventionally defined in terms of five parameters but it actually has only four independent 
parameters. Furthermore, the Wei |2] and Tietz potential functions [3] are exactly same 
solvable empirical functions. 
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Wang et al [4] also generated improved expressions for two versions of the Schioberg po¬ 
tential energy function which are the Rosen-Morse and Manning-Rosen potential functions. 
By choosing the experimental values of the dissociation energy, equilibrium bond length 
and equilibrium harmonic vibrational frequency as inputs, the authors obtained the average 
deviations of the energies calculated with the potential model from the experimental data 
for five diatomic molecules, and find that no one of six three-parameter empirical poten¬ 
tial energy functions is superior to the other potentials in fitting experimental data for all 
molecules examined. 

All these efforts were made in an attempt to find a most suitable molecular potential 
in its description of diatomic molecules. Following Refs. Bffl, we suggest sTW as a 
modification for the TW [2-8]. This potential can be written as 

- l) e - 6 W~re) _ ( c 2 _ X ) e -2 b h {r-r e ) 

(1 - c h e- b ^ r ~ r ^) 2 



V{r) = V e 


2 (c h 


where bh = 1 — Ch ), r e is the molecular bond length, (3 is the Morse constant, V e is 

the potential well depth and Ch is an optimization parameter obtained from ab initio or 
Rydberg-Klein-Rees (RKR) intramolecular potentials, r is the internuclear distance. When 
the potential constant approaches zero, i.e. Ch —> 0, the sTW potential reduces to the Morse 
potential [9]. This potential is just the TW potential shifted by dissociation energy D e . The 
shape of this potential is shown in figure [Ur, for different molecules. 

Figure [ljo compare between TW diatomic molecular potential, sTW diatomic molec¬ 
ular potential and the Morse potential using the parameters set for H 2 (A^E+J diatomic 
molecule. As it can be seen from this plot, the shifted Tietz-Wei and the Morse potentials 
are very close to each other for large values of r in the regions r « r e and r > r e , but they 
are very different at r « 0. This implies that the shifted Tietz-Wei potential is as good 
as traditional Morse potential and better than the Tietz-Wei potential in stimulating the 
atomic interaction for diatomic molecules. 

The scheme of our presentation is as follows. In the next section we give basic ingredient 
of exact quantization rule and all necessary formulas for our calculations. We solve the 
radial Schrodinger equation for the sTW and also obtain the rotational-vibrational energy 
spectrum for some diatomic molecules in section 3. Finally, results and conclusions are 
presented in section 4. 
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Figure 1: (a) Shape of Tietz-Wei diatomic molecular potential for different diatomic molecules, (b) Shape 
of the sTW, TW and Morse oscillator potentials for H 2 (X 1 E+) diatomic molecule with Ch = 0.170066, 
b h = 1.61890l-\ r e = 0.7416H and V e = 38318cm' 1 . 

2 A Brief Review of the Exact Quantization Rule 


Here we give a brief review on the EQR. The details can be found in refs. [10-24], In 
2005, Ma and Xu [I % IB] by carefully studying one-dimensional Schrodinger equation, have 
extended results to three-dimensional case by simply making the replacements x —> r and 
V(x) K/,(r): 


m 


k(r)dr = Nn + 


rn 

/ 0 ( r ) 

dk(r ) 

d(j)(r) 

L n } 

dr 

dr 


dr, k(r)= \Jj^l E ~V eff (r)}, (2) 


where and xb are two turning points determined by E = V e ff(r). The N = n + 1 is 
the number of the nodes of </>(r) in the region E n t = V e ff(r) and is larger by one than 
the number n of the nodes of wave function ij)(r). The first term Nn is the contribution 
from the nodes of the logarithmic derivative of wave function, and the second is called the 
quantum correction. 

In this approach, the energy spectrum equation is obtained by solving the two integrals 
involved in equation (J2]). This quantization rule has been used in many physical systems 
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to obtain the exact solutions of many exactly solvable quantum systems [10-24], EQR is a 
very important foundation to proper quantization rule (PQR) [25] . 


3 The Energy Spectrum 


To study any quantum physical model characterized by the diatomic molecular potential 
given by equation (JT]) , we need to solve the following Schrodinger equation for spherically 
symmetric potential in any arbitrary dimensional space: 
d 2 D - 1 d £{£ + D - 2) 2/i 


(■ E n i - V(r)) ipn,£, m (r : Q D ) = 0 . 


( 3 ) 


d 2 R n£ (r) 

2/i 

p . V 

"2 (c h - l)e~ bh d~ r «) - (, c 2 - l) e -2 b h {r-r e y 

2 d 2 -i)R 2 ] 

dr 2 

h 2 

y e 

1 - c h e~ bh d- r <d 

2/xr 2 


. dr 2 r dr r 2 % 

Now, by defining the wavefunction , Dd ) as r^ D ^ 2 R n£ (r)Y £m (6 , 0) and taking V(r) 

as the sTW diatomic molecular potential, the radial part of equation ([3]) can be found as 

fsfl _ If * 2 1 

Rni(r) = 0 , 

( 4 ) 

where n, £ and E n£ denote the principal quantum numbers, orbital angular momentum 
numbers and the bound state energy eigenvalues of the system under consideration (i.e., 

Eni < 0 ), respectively. The parameter 77 = £ + ^(D — 2) which is a linear combination of 
the spatial dimensions D and the angular momentum quantum number £. It is well known 
that for £ = 0, problem in the form (j3J) is exactly solvable. But for £ 7 ^ 0, it is not due 
to the centrifugal barrier. Therefore, in order to solve the above equation for £ 7 ^ 0 states, 
Hamzavi et al | 6 ] found that the following formula is a good approximation scheme to deal 
with the centrifugal barrier: 


9 1 Dq + D\- - , 7 

i — CfiC~° h ^ r — Te ' 


0 -b h (r-r e ) 


+ D-2 


?-2b h (r-r e ) 


(1 - c h e~ bh d- r -)y 


( 5 ) 


with 


D 0 = 1 - -(1 - c h )(3 + c h ) + -^(1 - c h ) 2 , lim D 0 = 1 - - + (6a) 

a or c h ^f 0 a cr 

D\ = -(1 - c h f{2 + c h ) - -^(1 - c h ) 3 , lim D 1 = — - (6b) 

a cr ch^-o a cr 

D ‘2 = --(1 - C h ) 3 (l + c h ) + -^(1 - c h ) A , lim D 2 = (6c) 

a cr c h -* 0 a cr 

Constant a = bh,r e has been introduced for mathematical simplicity. Now, by inserting this 

approximation into equation ([4]) and then introducing a new transformation of the form 
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r —* ( = L y- S - through the mapping function ( = /(r) with r in the domain [0, oo) or ( in 
the domain [—1, oo], we obtain the following second order differential equation: 


+ % Pm - V ef/ (()1 JWC) = 0. 


rl d( 

Ve„( C) = 


h z 
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e(e+i)D 2 

r-2 


with 

+ %£(c h - l) 2 




C/, 




Ch) 


T (e“ c 


.A 

J 7 



ZV^ 

2/i 

irf 

2/i 



+ R e (c„ - l) 2 . 


D x n 2 

2/i 


+ 214(0/, — 1), 



(7) 


If we define <f = e „^_ Cfa , then we can obtain the two turning points 4 and 4 as well as their 
sum and product properties by solving 14//(C) — E n i — 0 or 14//0?) — = 0 as: 


4 = - 7 ^ - - 4F(4l - and 4 = £= + V& ~ 4F(4l - E n ,\ 8) 


2J 7 2J 7 

B 


27 27 


^>a H - Q> 1 


A — E< 








and %) = \J^r h (5 - 4) (£ - 4)] 1/2 


Now, we can write the non-linear Riccati equation for the ground state is as 

- Q dl±£4; (f) + ^ (f) + ?§ [ Sm - K//(f)] = 0 (9) 

1 e fl 

Since the logarithmic derivative 4o(?) for the ground state has one zero and no pole, it has 
to take the linear form in 4 Thus, we assume the following solution for the ground states 


0o(?) — 4 + Bq 


( 10 ) 


By putting equation (TT0|) into (J9]), we can solve the non-linear Riccati equation. After proper 
comparison, it is straightforward to obtain the ground state energy and values of A and B 
as 


h 2 A 2 

E 0 £ = A -with A 
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pi B-T/cj , B 
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Since we now have all basic ingredient required to perform our calculations, thus, we proceed 
to calculating integrals (J2J) 

_ _re r ob 1 2 fiT [- (q - q a ) (q - 4)] 1/2 

' V - 1 - 1 ! 01 J^Oa ' 

7T7V /2/iJ 7 
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a 
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a 
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( 12 ) 
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where we have used the following standard integral 


" Xb y/(:x a -x){x- x b ) 
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J Xa ^ + QZ) 

Furthermore, we can hnd the integral at the right hand side as 


(13) 
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Using the results in equations (fl4jl and dT3|) with equation (J2J), we can hnd the energy energy 
eigenvalues equation for the sTW diatomic molecular potential as 

(5 + n ) 2 + ^ q2c p {Di Ch ~ Jh) + (l - jr 


(14) 
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In three-dimensions (D — 3), it can be reduced to the form 


2*2 




h £{£ + 1) A) a 2 ti 


2/ir, 
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(16) 


1 1 


with 5 = - + -\l 1 + — 


4 The Eigenfunctions 


4 /ZV^ + l) 2/rKr e 2 


cr 


+ 


a 2 r 




For the sake of completeness, we study the corresponding wavefunctions for this potential. 
For this purpose we introduce a new transformation of the form t = e~ bh ^ r ~ r ^ E (e a ,0) 
in equation (17|h which maintained the finiteness of the transformed wave functions on the 
boundary conditions. Thus, we can find 


d 2 U ni (t ) 1 dU nl (t) 


dt 2 
-2 c h 


t dt ' t 2 (1 — Cht ) 2 
2 iir 2 e (E nl + V e ) (r) 2 -D 


Dn 


2nr 2 e E n i _ (■ fj 2 - |) 
h 2 a 2 a 2 

4/rr 2 K (fj 2 - |) 


a z 


+ 
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a 2 h 2 

2 ttr 2 e (E ne + V e ) (v 2 - ±) 
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a 2 h 2 
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h 2 a 2 

(fi 1 - \) 


a z 


D i 


cr 


( d ^- d 2 )- 2 -0 

n a 2 


(17) 


t 2 > U„i(t) = 0. 


Following the procedure described in ref. [26], we can write the solution U n e(t) in terms of 
hypergeometric polynomials and thus, the wave function takes the form 


U n i(Q = N n ie r ° c (l-c h e “ c )' 1 2*1 (- 11 ,n + 2(p + 4 ); 2p + l,c h e a< ) , 


(18) 


with 


V = 


9 0 J-'nl 


a z 


k 2 a 2 


and? = ^ 1 + \/ 1 + 5 L 


V 2 ~\\ n , 2/rr 2 K 


cr 


D 2 + 


h 2 a 2 


(1 - c h ) 2 
(19) 


where N n i is the normalization constant. For a further detail on the calculation of similar 
potential models solved using formula method, one is advised to refer to other work 


5 Results and Conclusions 


In this study, in an attempt to find a more suitable potential that stimulate the atomic 
interaction in diatomic molecules, we suggested sTW diatomic molecular potential as a 
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Table 1: Model parameters of the diatomic molecules studied in the present work. 


Molecules (states) 

Ch 

m/io 23 (g) 

b^A- 1 ) 

r e (A) 

D(cm - 1 ) 

NO (a 4 IIi) 

0.0082003 

1.249 

2.408413 

1.451 

16361 

NO (b 2 n r ) 

-0.482743 

1.249 

3.42650 

1.428 

22722 

NO (L ,2 0) 

-0.073021 

1.249 

2.73796 

1.451 

14501 

NO (6 4 S“) 

-0.085078 

1.249 

3.01538 

1.318 

21183 

ici (xis+) 

-0.086212 

4.55237 

2.008578 

2.3209 

17557 

IC1 (A 3 IIi) 

-0.167208 

4.55237 

2.542557 

2.6850 

3814.7 

ici (A ,3 n 2 ) 

-0.157361 

4.55237 

2.373450 

2.6650 

4875 


Table 2: The bound states energy eigenvalues (D = 2 and 3) for set of diatomic molecules 
for various n and rotational i quantum numbers in sDF diatomic molecular potential. 


D 

n 

l 

NO (a 4 nd 

NO (B 2 n r ) 

NO (L' 2 <t>) 

NO (6 4 £-) 

ICI(X 4 S+) 

ici M 3 ni) 

ICI (A l3 n 2 ) 


0 

0 

-1.971298585 

-2.88233770 

-1.855428285 

-2.695026855 

-2.200695845 

-0.4861418795 

-0.618432865 


1 

0 

-1.859308585 

-3.01386965 

-1.972859085 

-2.834614855 

-2.248819245 

-0.5129225795 

-0.646829015 



1 

-1.859148585 

-3.01373250 

-1.972723585 

-2.834451355 

-2.248804845 

-0.5129118495 

-0.646818105 


2 

0 

-1.750608585 

-3.14759050 

-2.093631485 

-2.977436555 

-2.297420945 

-0.5403157895 

-0.675774025 



1 

-1.750468585 

-3.14745280 

-2.093493685 

-2.977270755 

-2.297406745 

-0.5403049695 

-0.675763045 



2 

-1.750098585 

-3.14703970 

-2.093081285 

-2.976773655 

-2.297363645 

-0.5402725095 

-0.675730195 

2 

3 

0 

-1.645178585 

-3.28347662 

-2.217724085 

-3.123472855 

-2.346500045 

-0.5683148895 

-0.705263165 



1 

-1.645068585 

-3.28333841 

-2.217584185 

-3.123304955 

-2.346485745 

-0.5683039695 

-0.705252125 



2 

-1.644678585 

-3.28292375 

-2.217165185 

-3.122801255 

-2.346442645 

-0.5682712195 

-0.705219005 



3 

-1.644058585 

-3.28223274 

-2.216466785 

-3.121961855 

-2.346370345 

-0.5682166495 

-0.705163785 


4 

0 

-1.543078585 

-3.42150558 

-2.345115985 

-3.272704955 

-2.396055245 

-0.5969133695 

-0.735291785 



1 

-1.542938585 

-3.42136686 

-2.344973985 

-3.272534755 

-2.396040945 

-0.5969023495 

-0.735280655 



2 

-1.542598585 

-3.42095075 

-2.344548485 

-3.272024355 

-2.395997445 

-0.5968693195 

-0.735247275 



3 

-1.541978585 

-3.42025729 

-2.343839285 

-3.271173455 

-2.395924945 

-0.5968142695 

-0.735191635 



4 

-1.541118585 

-3.41928647 

-2.342846285 

-3.269982555 

-2.395823745 

-0.5967371995 

-0.735113745 


0 

0 

-1.971278585 

-2.88230356 

-1.855395085 

-2.694986655 

-2.200692145 

-0.4861392095 

-0.618430175 


1 

0 

-1.859278585 

-3.01383536 

-1.972825385 

-2.834574055 

-2.248815545 

-0.5129198795 

-0.646826285 



1 

-1.859008585 

-3.01356106 

-1.972554585 

-2.834247155 

-2.248787045 

-0.5128984495 

-0.646804535 


2 

0 

-1.750568585 

-3.14755607 

-2.093596885 

-2.977395155 

-2.297417445 

-0.5403130995 

-0.675771295 



1 

-1.750288585 

-3.14728068 

-2.093321985 

-2.977063955 

-2.297388745 

-0.5402914395 

-0.675749345 



2 

-1.749808585 

-3.14672991 

-2.092771985 

-2.976400855 

-2.297331445 

-0.5402481695 

-0.675705525 

3 

3 

0 

-1.645158585 

-3.28344206 

-2.217689085 

-3.123430855 

-2.346496645 

-0.5683121595 

-0.705260395 



1 

-1.644908585 

-3.28316563 

-2.217409585 

-3.123095155 

-2.346467645 

-0.5682903295 

-0.705238345 



2 

-1.644418585 

-3.28261279 

-2.216850885 

-3.122423255 

-2.346410245 

-0.5682466595 

-0.705194145 



3 

-1.643678585 

-3.28178358 

-2.216012785 

-3.121415755 

-2.346323545 

-0.5681811695 

-0.705127915 


4 

0 

-1.543028585 

-3.42147091 

-2.345080285 

-3.272662455 

-2.396051945 

-0.5969106095 

-0.735288995 



1 

-1.542798585 

-3.42119348 

-2.344796885 

-3.272322355 

-2.396022645 

-0.5968885995 

-0.735266735 



2 

-1.542328585 

-3.42063869 

-2.344229385 

-3.271641355 

-2.395964745 

-0.5968445495 

-0.735222225 



3 

-1.541588585 

-3.41980654 

-2.343378185 

-3.270620455 

-2.395877945 

-0.5967784695 

-0.735155455 



4 

-1.540608585 

-3.41869709 

-2.342243285 

-3.269259655 

-2.395761945 

-0.5966903895 

-0.735066435 














Table 3: The bound states energy eigenvalues D = 4 and 5) for set of diatomic molecules 
for various n and rotational i quantum numbers in sDF diatomic molecular potential. 


D 

n 

t 

NO (a 4 nd 

NO (B 2 n r ) 

NO (L' 2 </>) 

NO (6 4 E") 

IC1(x 4 s+) 

IC1 (zi 3 ni) 

ici h4' 3 n 2 ) 


0 

0 

-1.971158585 

-2.88220112 

-1.855295085 

-2.694865755 

-2.200681645 

-0.4861312495 

-0.618422105 


1 

0 

-1.859148585 

-3.01373250 

-1.972723585 

-2.834451355 

-2.248804845 

-0.5129118495 

-0.646818105 



1 

-1.858778585 

-3.01332103 

-1.972317885 

-2.833961155 

-2.248762245 

-0.5128796695 

-0.646785525 


2 

0 

-1.750468585 

-3.14745280 

-2.093493685 

-2.977270755 

-2.297406745 

-0.5403049695 

-0.675763045 



1 

-1.750098585 

-3.14703970 

-2.093081285 

-2.976773655 

-2.297363645 

-0.5402725095 

-0.675730195 



2 

-1.749438585 

-3.14635128 

-2.092393885 

-2.975945655 

-2.297292045 

-0.5402183995 

-0.675675405 

3 

3 

0 

-1.645068585 

-3.28333841 

-2.217584185 

-3.123304955 

-2.346485745 

-0.5683039695 

-0.705252125 



1 

-1.644678585 

-3.28292375 

-2.217165185 

-3.122801255 

-2.346442645 

-0.5682712195 

-0.705219005 



2 

-1.644058585 

-3.28223274 

-2.216466785 

-3.121961855 

-2.346370345 

-0.5682166495 

-0.705163785 



3 

-1.643188585 

-3.28126535 

-2.215488785 

-3.120786155 

-2.346269545 

-0.5681402395 

-0.705086495 


4 

0 

-1.542938585 

-3.42136686 

-2.344973985 

-3.272534755 

-2.396040945 

-0.5969023495 

-0.735280655 



1 

-1.542598585 

-3.42095075 

-2.344548485 

-3.272024355 

-2.395997445 

-0.5968693195 

-0.735247275 



2 

-1.541978585 

-3.42025729 

-2.343839285 

-3.271173455 

-2.395924945 

-0.5968142695 

-0.735191635 



3 

-1.541118585 

-3.41928647 

-2.342846285 

-3.269982555 

-2.395823745 

-0.5967371995 

-0.735113745 



4 

-1.540048585 

-3.41803836 

-2.341569785 

-3.268451355 

-2.395693345 

-0.5966380895 

-0.735013565 


0 

0 

-1.970998585 

-2.88203041 

-1.855128785 

-2.694664455 

-2.200663845 

-0.4861179795 

-0.618408615 


1 

0 

-1.859008585 

-3.01356106 

-1.972554585 

-2.834247155 

-2.248787045 

-0.5128984495 

-0.646804535 



1 

-1.858508585 

-3.01301246 

-1.972013385 

-2.833593555 

-2.248729945 

-0.5128555495 

-0.646761065 


2 

0 

-1.750288585 

-3.14728068 

-2.093321985 

-2.977063955 

-2.297388745 

-0.5402914395 

-0.675749345 



1 

-1.749808585 

-3.14672991 

-2.092771985 

-2.976400855 

-2.297331445 

-0.5402481695 

-0.675705525 



2 

-1.749048585 

-3.14590380 

-2.091946885 

-2.975407055 

-2.297245545 

-0.5401832395 

-0.675639815 

4 

3 

0 

-1.644908585 

-3.28316563 

-2.217409585 

-3.123095155 

-2.346467645 

-0.5682903295 

-0.705238345 



1 

-1.644418585 

-3.28261279 

-2.216850885 

-3.122423255 

-2.346410245 

-0.5682466595 

-0.705194145 



2 

-1.643678585 

-3.28178358 

-2.216012785 

-3.121415755 

-2.346323545 

-0.5681811695 

-0.705127915 



3 

-1.642668585 

-3.28067802 

-2.214895385 

-3.120072855 

-2.346208245 

-0.5680938395 

-0.705039565 


4 

0 

-1.542798585 

-3.42119348 

-2.344796885 

-3.272322355 

-2.396022645 

-0.5968885995 

-0.735266735 



1 

-1.542328585 

-3.42063869 

-2.344229385 

-3.271641355 

-2.395964745 

-0.5968445495 

-0.735222225 



2 

-1.541588585 

-3.41980654 

-2.343378185 

-3.270620455 

-2.395877945 

-0.5967784695 

-0.735155455 



3 

-1.540608585 

-3.41869709 

-2.342243285 

-3.269259655 

-2.395761945 

-0.5966903895 

-0.735066435 



4 

-1.539388585 

-3.41731030 

-2.340824985 

-3.267557955 

-2.395617345 

-0.5965802895 

-0.734955135 


modification for the TW diatomic molecular potential. The bound state solution of this 
potential has been found in an arbitrary D-dimension via the improved exact quantization 
rule. 

Further, using the spectroscopic parameters presented in table [T| which are taken from 
ref. |8j, we computed rotational-vibrational energy spectrum of some diatomic molecules 
in 2,3,4,5-dimensions. The results are presented in tables [2] and [3j In our numerical com¬ 
putations, we have used the following conversions: lamu = 931.494028MeI//c 2 , 1cm -1 = 
1.239841875 x 10~ 4 eV, and he = 1973.29 eVA. 
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From equation (1T5]) . it can be seen that two interdimensional states are degenerate 
whenever ( n,£,D ) —>■ (n, £ ± 1,D =F 2) E^ ( = E^j^f . Thus, a knowledge of E^ e for 

D = 2 to 5 provides the information necessary to find E® e for other higher dimensions. 
For example, Eq 2 ] = EqI = . This is the same transformational invariance 

described for bound states of free atoms and molecules [2711281129] and demonstrates the 
existence of interdimensional degeneracies among states of the confined Hulthen potential. 

The advantage of the approach employed in this study is that it gives the eigenvalues 
through the calculation of two integral given by equation 0 and solving the resulting 
algebraic equation. Firstly, we can easily obtain the quantum correction by only considering 
the solution of the ground state of the quantum system since it is independent of the number 
of nodes of the wave function for exactly solvable quantum system. The general expressions 
obtained for the energy eigenvalues and wave functions can be easily reduced to the 3D 
space (D = 3) and for s-wave (i.e. £ = 0 state). The EQR produce as good results as the 
PQR, however the procedure followed using PQR is more shorter and quick. 
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